Cardiac anisotropy in boundary-element models for the electrocardiogram

The boundary-element method (BEM) is widely used for electrocardiogram (ECG) simulation. Its major disadvantage is its perceived inability to deal with the anisotropic electric conductivity of the myocardial interstitium, which led researchers to represent only intracellular anisotropy or neglect anisotropy altogether. We computed ECGs with a BEM model based on dipole sources that accounted for a “compound” anisotropy ratio. The ECGs were compared with those computed by a finite-difference model, in which intracellular and interstitial anisotropy could be represented without compromise. For a given set of conductivities, we always found a compound anisotropy value that led to acceptable differences between BEM and finite-difference results. In contrast, a fully isotropic model produced unacceptably large differences. A model that accounted only for intracellular anisotropy showed intermediate performance. We conclude that using a compound anisotropy ratio allows BEM-based ECG models to more accurately represent both anisotropies

3D finite element electrical model of larval zebrafish ECG signals

Assessment of heart function in zebrafish larvae using electrocardiography (ECG) is a potentially useful tool in developing cardiac treatments and the assessment of drug therapies. In order to better understand how a measured ECG waveform is related to the structure of the heart, its position within the larva and the position of the electrodes, a 3D model of a 3 days post fertilisation (dpf) larval zebrafish was developed to simulate cardiac electrical activity and investigate the voltage distribution throughout the body. The geometry consisted of two main components; the zebrafish body was modelled as a homogeneous volume, while the heart was split into five distinct regions (sinoatrial region, atrial wall, atrioventricular band, ventricular wall and heart chambers). Similarly, the electrical model consisted of two parts with the body described by Laplace’s equation and the heart using a bidomain ionic model based upon the Fitzhugh-Nagumo equations. Each region of the heart was differentiated by action potential (AP) parameters and activation wave conduction velocities, which were fitted and scaled based on previously published experimental results. ECG measurements in vivo at different electrode recording positions were then compared to the model results. The model was able to simulate action potentials, wave propagation and all the major features (P wave, R wave, T wave) of the ECG, as well as polarity of the peaks observed at each position. This model was based upon our current understanding of the structure of the normal zebrafish larval heart. Further development would enable us to incorporate features associated with the diseased heart and hence assist in the interpretation of larval zebrafish ECGs in these conditions

Impact of Sarcoplasmic Reticulum Calcium Release on Calcium Dynamics and Action Potential Morphology in Human Atrial Myocytes: A Computational Study

Electrophysiological studies of the human heart face the fundamental challenge that experimental data can be acquired only from patients with underlying heart disease. Regarding human atria, there exist sizable gaps in the understanding of the functional role of cellular Ca2+ dynamics, which differ crucially from that of ventricular cells, in the modulation of excitation-contraction coupling. Accordingly, the objective of this study was to develop a mathematical model of the human atrial myocyte that, in addition to the sarcolemmal (SL) ion currents, accounts for the heterogeneity of intracellular Ca2+ dynamics emerging from a structurally detailed sarcoplasmic reticulum (SR). Based on the simulation results, our model convincingly reproduces the principal characteristics of Ca2+ dynamics: 1) the biphasic increment during the upstroke of the Ca2+ transient resulting from the delay between the peripheral and central SR Ca2+ release, and 2) the relative contribution of SL Ca2+ current and SR Ca2+ release to the Ca2+ transient. In line with experimental findings, the model also replicates the strong impact of intracellular Ca2+ dynamics on the shape of the action potential. The simulation results suggest that the peripheral SR Ca2+ release sites define the interface between Ca2+ and AP, whereas the central release sites are important for the fire-diffuse-fire propagation of Ca2+ diffusion. Furthermore, our analysis predicts that the modulation of the action potential duration due to increasing heart rate is largely mediated by changes in the intracellular Na+ concentration. Finally, the results indicate that the SR Ca2+ release is a strong modulator of AP duration and, consequently, myocyte refractoriness/excitability. We conclude that the developed model is robust and reproduces many fundamental aspects of the tight coupling between SL ion currents and intracellular Ca2+ signaling. Thus, the model provides a useful framework for future studies of excitation-contraction coupling in human atrial myocytes

Mathematical Modeling and Simulation of Ventricular Activation Sequences: Implications for Cardiac Resynchronization Therapy

Next to clinical and experimental research, mathematical modeling plays a crucial role in medicine. Biomedical research takes place on many different levels, from molecules to the whole organism. Due to the complexity of biological systems, the interactions between components are often difficult or impossible to understand without the help of mathematical models. Mathematical models of cardiac electrophysiology have made a tremendous progress since the first numerical ECG simulations in the 1960s. This paper briefly reviews the development of this field and discusses some example cases where models have helped us forward, emphasizing applications that are relevant for the study of heart failure and cardiac resynchronization therapy

25th Annual Computational Neuroscience Meeting: CNS-2016

Abstracts of the 25th Annual Computational Neuroscience Meeting: CNS-2016 Seogwipo City, Jeju-do, South Korea. 2–7 July 201

Scalable Heterogeneous CPU-GPU Computations for Unstructured Tetrahedral Meshes

A recent trend in modern high-performance computing environments is the introduction of powerful, energy-efficient hardware accelerators such as GPUs and Xeon Phi coprocessors. These specialized computing devices coexist with CPUs and are optimized for highly parallel applications. In regular computing-intensive applications with predictable data access patterns, these devices often far outperform CPUs and thus relegate the latter to pure control functions instead of computations. For irregular applications, however, the performance gap can be much smaller and is sometimes even reversed. Thus, maximizing the overall performance on heterogeneous systems requires making full use of all available computational resources, including both accelerators and CPUs

Scalable Heterogeneous CPU-GPU Computations for Unstructured Tetrahedral Meshes

A recent trend in modern high-performance computing environments is the introduction of powerful, energy-efficient hardware accelerators such as GPUs and Xeon Phi coprocessors. These specialized computing devices coexist with CPUs and are optimized for highly parallel applications. In regular computing-intensive applications with predictable data access patterns, these devices often far outperform CPUs and thus relegate the latter to pure control functions instead of computations. For irregular applications, however, the performance gap can be much smaller and is sometimes even reversed. Thus, maximizing the overall performance on heterogeneous systems requires making full use of all available computational resources, including both accelerators and CPUs

Heart Closed-Loop Model for the Assessment of Cardiac Pacing

We developed a closed-loop model of cardiac stimulation using a finite element model of the whole-heart embedded in the torso that is proposed as an useful tool for pacemaker design and testing. The electrical activity of the cardiac tissue is reproduced with a bidomain model incorporated with the FitzHugh-Nagumo equations. The finite element model is developed in Comsol Multiphysics, both in two and in three dimensions and then exported in Simulink environment where the pacemaker algorithm is implemented. To validate the model, we chose a demand inhibited pacemaker, which stimulates the myocardium only if the intrinsic activity of the heart is not revealed, but every type of pacemaker can be simulated. The model generates a controlled spontaneous activation in the sinoatrial node and it is also able to reproduce realistic electrocardiographic signals and the effects that the stimulation has on them